Multiple Glasses in Asymmetric Binary Hard Spheres 
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Multiple distinct glass states occur in binary hard-sphere mixtures with constituents of very dis- 
parate sizes according to the mode-coupling theory of the glass transition (MCT), distinguished by 
considering whether small particles remain mobile or not, and whether small particles contribute 
significantly to perturb the big-particle structure or not. In the idealized glass, the four different 
glasses are separated by sharp transitions that give rise to higher-order transition phenomena in- 
volving logarithmic decay laws, and to anomalous power-law-like diffusion. The phenomena are 
argued to be expected generally in glass-forming mixtures. 

PACS numbers: pacs tbd 



Many glass formers and virtually all simple model sys- 
tems for slow dynamics are mixtures of some sort. Close 
to a glass transition generic mixing effects appear, such 
as dramatic changes in viscosity induced by small com- 
position changes [J, Q , that are relevant for applications 
and may help to shed light on the microscopic processes 
driving glass formation. 

Even more interesting is the possibility to form qual- 
itatively distinct types of glass, depending on mixture 
composition and constituents. Taking for example bi- 
nary mixtures with sufficiently disparate constituents, a 
glass can form where some (slow) species freeze, but a fast 
component is able to diffuse through the voids left in the 
amorphous packing. This scenario is particularly relevant 
for transport through heterogeneous disordered media 0- 
7] or glassy ion conductors [1, 0. The simplest model 
are binary hard-sphere mixtures with large size dispar- 
ity, where experiments on colloidal suspensions indeed 
found, depending on relative concentration, a partially 
frozen "single glass" with mobile small particles, sepa- 
rate from a "double glass" where both particle species 



freeze 13, 11 



In mixtures of star polymers [12H14] yet 
another kind of glass emerged, termed "asymmetric" be- 
cause it is characterized by few big particles frozen in a 
small-particle matrix, rendering the big-particle nearest- 
neighbor cages highly nonspherical. It was, however, ar- 
gued to be a hallmark of the ultra-soft interactions typi- 
cal for the star polymers. 

Another kind of glass intuitively argued for is the 
"attractive glass" famous from colloid-polymer mixtures 
where free polymer induces depletion attraction among 
the colloids. If that attraction is weak, the glass that 
forms is essentially hard-sphere like or "repulsive" , while 
at sufficiently strong and short-ranged attraction, a new 
glass driven by bonding and not nearest-neighbor cage- 
ing appears. Based on extensively tested predictions of 
the mode-coupling theory of the glass transition (MCT) 
for a square- well model system la- 13], one expects the 
two glasses to be separated by a glass-glass transition 



crossing which, for example, the elastic moduli of the 
glass exhibit sharp changes. Considering the generality 
of the depletion-interaction mechanism 2l| , one may in- 
deed expect a similar glass-glass transition to be present 
in binary mixtures quite generically (lest the relevant pa- 
rameter space cannot be explored). As such transitions 
typically involve endpoint singularities that give rise to 
universal logarithmic decay laws for the time-dependent 
correlation functions [22| , one anticipates regions of mix- 
ture composition where these peculiar decay laws can be 
found. Indeed, they have been reported in computer sim- 
ulation of soft-sphere mixtures [231 ] . 

Here I demonstrate, that already the simplest glass- 
forming mixture model, the binary hard-sphere mixture, 
allows to identify four qualitatively different glassy states 
separated by well-defined transitions. The transition di- 
agram lends itself to an intuitive classification: one has 
to distinguish (i) whether small particles remain mobile 
in the glass or not, and (ii) whether the structure of the 
big-particle glass is or is not significantly disturbed by 
the small particles. Combining these two choices each 
gives four possible types of glass. 

Calculations are based on MCT supplemented by the 
Percus-Yevick (PY) approximation for the static struc- 
ture factor [2J-|26[, but following the physically plausi- 
ble classification, the results can be expected to hold 
rather generally for binary mixtures whose parameters 
can be tuned widely enough. The possible interplay with 
equilibrium phases is ignored here. The resulting glass- 
transition diagram is a unique prediction of MCT, dis- 
tinct from other theories that have been put forward [13] , 
and testable in simulation or experiment. 

Numerical calculations follow the method of Ref. (2f| . 
MCT takes partial static structure factors S a p(q) = 
(Qa( ( l)Ql3{ < l)) as input (Greek indices label species, and 
static triplet correlations are neglected here). Here (•) 
is the canonical average, and Q a {<t) = Y2k ex v[i < Fk,a] 
is the number-density fluctuation of species a (the sum 
runs over all particle positions r of that type). MCT 
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x (vol-% small particles) 

FIG. 1. (a) Glass transition diagram for a model binary hard- 
sphere mixture obtained from mode-coupling theory. Solid 
lines show the glass-transition packing fraction tp c (S, x) as a 
function of small-particle volume concentration, x, for various 
size ratios 8 as indicated. The dashed line indicates a separate 
localization transition tp a c (8,x) for the small species for 8 — 
0.2. (b) MCT exponent parameter \(x) for cuts of constant 
8. Solid and open symbols refer to the small-! and large- a: 
branches, respectively. 

predicts collective dynamical density correlation func- 
tions 3> a p(q,t) — (Qa(q,t)gp(q)), where the time evo- 
lution is given by the implicit dependence of g a (q) 
on the trajectories. The long-time limit of that ma- 
trix, F(q) — lim^oo t), distinguishes ergodic liq- 
uid states, F(q) = 0, from nonergodic ones, where the 
correlation matrix decays not to zero but to a positive 
definite matrix, F(q) >- 0. Such states are identified 
as (idealized) glasses within MCT. Standard MCT glass 
transitions are characterized as bifurcations of the A2 
type where F(q) jumps discontinuously. Note that such 
jumps can also occur inside the glass. The F(q) can 
be found numerically [28| on a finite wave-number grid, 
q l = (i+l/2)Aq with i = 1, . . . N = 500 and Ag = 0.4 for 
most calculations here. Small size ratios require higher 
large-g cutoff, N — 1000 was hence used for S < 0.2. 

Results are reported as cuts through the parameter 
space at constant size ratio 5 = d s /d\ < 1, where d a are 
the hard-sphere diameters (a = s, 1 for small and large). 
The two remaining parameters then are the number den- 
sities p a , conveniently expressed as a total packing frac- 



tion of the spheres, ip = J2 a fa = Or/6) Y, a Pad%, and 
a concentration (by volume) of small spheres, i: — tp s /tp- 
To top panel of Fig. Q] illustrates the glass-transition di- 
agram for size ratios 0.2 < S < 0.5; larger 8 are topo- 
logically the same as 8 = 0.5 (25|. For all compositions 
at fixed 8 > 0.42, the binary-mixture glass is separated 
from the liquid by a smooth line of ordinary A2 transi- 
tions. In contrast, the curve for, say, 8 = 0.35 demon- 
strates the emergence of a glass-glass transition: the 
liquid-glass transition splits into two lines that no longer 
join smoothly, but cross discontinuously at some point, 
at which one of the lines stops. The other continues in- 
side the glass until it terminates (the corresponding jump 
in F(q) disappears) at a higher-order transition point of 
type ^3 . The glass-glass transition signals that there are 
two competing arrest mechanisms at work: cageing can 
either be dominated by large particles or by small ones. 
As soon as the ratio of relevant length scales (essentially 
5) becomes sufficiently distinct from unity, the two mech- 
anisms give rise to glasses with differing intrinsic scales 
(the localization length, or typical cage size) such that a 
sharp distinction becomes possible. As pointed out, one 
way to distinguish these glasses experimentally would be 
a marked difference in elastic properties. The structure 
of the large- x glass bears resemblance to certain types of 
sweets [27[ such as Italian torroncino. The distinction is 
only strict for size ratios 8 < 5 C . At S c , the glass-glass 
transition endpoint coincides with the crossing, resulting 
in an A4 singularity. Mathematically, these higher-order 
singularities are identical to those in one-component de- 
scriptions of colloid-polymer mixtures (22j , and the same 
kind of asymptotic expansions predicting the the appear- 
ance of logarithmic decay laws and their precursors apply 
here, too. In the present calculations, S c « 0.41. Since 
the states x = and x — 1 are identical up to a rescal- 
ing of lengths with 6, the .F-versus-g curve is broader for 
x — > 1 than for x — > 0, and generically higher at fixed q. 
Hence the large- a; portion of the transition extends into 
and demarks the mechanically stiffcr glass. 

The position of higher-order singularities is revealed by 
looking at MCT's exponent parameter A. To each tran- 
sition point, a particular value of 0.5 < A < f is assigned 
in the theory, which determines the power-law exponents 
valid for the asymptotic expansion of correlation func- 
tions close to the transition. At A = 1, these power laws 
cease to be valid and are replaced by logarithmic laws. 
As seen in the lower panel of Fig. [TJ the values of A in- 
deed rise to unity when approaching the endpoints of the 
large- a; glass-transition lines. At 8 > S c , A still rises as 
a precursor corresponding to a notable decrease in the 
asymptotic exponents; at 8 — 5 C , the A(£)-curve reaches 
unity for the first time before splitting into two branches. 

MCT predicts that at all the At transition points 
discussed so far, all collective partial-density correla- 
tion functions show a simultaneous jump in their long- 
time limits and hence all are nonergodic [281 ] . However, 
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this does not need to hold for the nonergodicity factors 
fail) = nm t^oo 0a(<ZiO associated with the self-motion 
of the species a, <j) s a (q 7 t) = (exp[iq(r s ^ a (t) - r SjQ (0))]). 
fs(g) can remain zero when crossing the glass transition 
line, if the small particles are below a certain size 8 8 c 
[U & A dashed line in Fig. Ufa) demarcs for 5 = 0.2 
the region where /* (q) — and hence the "single glass" . 
The collective small-small density correlation function 
remains nonergodic also in the single glass, in distinct 
disagreement with the SCGLE theory of Rcf. 
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The nature of the transition lines shown in dashed in 
Fig. [1] is notably different from the At glass transition 
lines shown as solid lines: they are so-called type A or 
localization transitions, where the nonergodicity factor 
fs(q) rises continually from zero upon crossing the transi- 
tion, as opposed to exhibiting a finite jump. This reflects 
the different nature of the two transition mechanisms: 
while ordinary glass transitions are driven by collective 
caging on local length scales, the localization transition 
contains a divergent two-point length scale (the localiza- 
tion length of the small particles) . This also means that 
the applicability of MCT to such transitions can be de- 
bated (even more than MCT is usually debated). The 
dashed line in Fig. QJa) is to be seen as a good qualita- 
tive approximation to the localization line [3l| , as it is ex- 
pected on physical grounds: increasing the density of the 
single glass, one will eventually reach a regime where the 
small particles either have too little void space to move, 
or are dense enough to form a glass on their own. For 
this reason, the single-glass region is bounded from below 
by the liquid-glass transition, from above by a localiza- 
tion transition, and bounded from the high-x side by the 
glass-glass transition in Fig. [TJa). This is in variance 
with Ref. [27[, where no glass-glass transition is found. 
The separate localization transition in our calculation oc- 
curs ony below some size ratio <5* strictly smaller than 
the 5 C where the higher-order singularity first occurs. 
The precursor to the localization transition is anoma- 
lous diffusion in the sense that the small-particle mean- 
squared displacement, 5r 2 {t) = ((f s , s (t) — rv s (0)) 2 }, ex- 
hibits pow er-law growth Sr 2 (t) oc t y with an exponent 
y < 1 [29(, instead of ordinary diffusion (y — 1), or the 
two-step glass-transition pattern comprising an interme- 
diate plateau (y ~ 0). This allows the existence of the 
"double-transition" scenario to be established, as done 
recently for simulations of a soft-sphere mixture 

So far, we have found three distinct glasses separated 
by well-defined transitions. The glass-glass transition 
is similar to the one found in colloid-polymer mixtures: 
both share the same mathematics, and both arise from 
a competition of arrest mechanisms on two sufficiently 
distinct length scales. Tempting as it may be [32j, the 
analogy is flawed, since in the colloid-polymer mixture 
the small component is always assumed to remain mobile, 
whereas at the glass-glass transition in Fig. [TJ the small 
component itself becomes the main glass former. 
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FIG. 2. Transition diagram as in Fig. [1] but for smaller S. 
Open, filled, and shaded symbols correspond to X(x) for the 
small-i, large-f, and intermediate-i branch of <j> c (x), respec- 
tively. Inset: localization length r\ for large and small parti- 
cles for (8,x) = (0.18,0.45) as a function of packing fraction. 



The analog of the attractive glass in the binary mixture 
has to be found inside the single glass region. Indeed, 
extending Fig. Q] to smaller 5, Fig. [2 a second glass- 
glass transition emerges: a second A4 point at roughly 
5* = 0.194 marks the onset of the higher-order singu- 
larity scenario for S < 5* at around x rs 0.3. Again, it 
indicates the discontinuous change of a big-particle domi- 
nated glass to one where the small particles set a relevant 
length scale, but now in the sense that the small parti- 
cles induce a strong depletion attraction while themselves 
remaining mobile. 

The different glasses are distinguished by the localiza- 
tion lengths r s c a of big and small particles - a measure of 
their cage size. The inset of Fig.^b) shows r s c a (ip) along 
a cut crossing all four transitions. As the repulsive single- 
glass is first entered, r s cl « 0.1, the Lindemann length for 
big particles. It discontinuously drops by about 5 enter- 
ing the attractive single-glass, as depletion forces move 
the large particles closer together. The small-particle 
r s c s remains infinite up to the localization transition, and 
shows signs of a continuous divergence there (where l 
is smooth). As the "torroncino" glass is entered, both 
r s c drop discontinuously again, and r s cs 0.1(5 indicates 



a small-particle glass. 
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To summarize, four different glasses are predicted to 
form in the simplest model of glass-forming mixtures, 
viz., binary hard-sphere mixtures, if both composition 
and size ratio are changed. The four glasses come in 
two categories, one ("double glass") where both species 
freeze simultaneously, and one ( "single glass" ) where the 
smaller component remains mobile inside the frozen en- 
vironment. Each of these comes in two variants, "repul- 
sive" and "attractive" , depending on whether the struc- 
ture is primarily driven by large-particle caging, or by 
the small-particle induced forces (be they arrested or 
not). The attractive double glass naturally explains low- 
big-particle coordination "asymmetric" cages thought to 
arise from ultra-soft potentials. Our calculations make 
it plausible to expect the existence of four glasses much 
more generically, as long as the relevant mixture param- 
eters cover a sufficiently wide range of states (which may, 
for example, also depend on dimensionality |26j). 

A unique prediction is that the four types of glass are 
separated by sharp transitions and accompanied by re- 
gions in parameter space where signatures of higher-order 
singularities or of anomalous diffusion may be seen. The 
latter was found in recent simulations [§] . Features of log- 
arithmic decay that likely are the signature of the first set 
of ^3 singularities have been reported [IH • The possibil- 
ity of attractive glasses emerging around S ~ 0.2 » 6* has 
been hinted at [33] . Further simulations and experiments 
are called for to conclusively test our predictions. 

Cautionary remarks should be added: first, the PY ap- 
proximation has known defects. Still it appears to pre- 
dict S(q) surprisingly sensibly The effects discussed 
above arise from an interplay of different length scales; 
this should be captured qualitatively correctly in PY. 
Second, experiment and simulation may require to resort 
to more complicated interaction potentials, to suppress 
equilibrium phase transitions that could otherwise inter- 
fere. Finally, MCT is driven to a regime where it may 
likely fail; our predictions thus pose a demanding test 
case for the theory. 
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